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Abstract 



Biological adhesion often involves several pairs of specific receptor-ligand 
molecules. Using rate equations, we study theoretically the rupture of such 
multiple parallel bonds under dynamic loading assisted by thermal activation. 
For a simple generic type of cooperativity, both the rupture time and force 
exhibit several different scaling regimes. The dependence of the rupture force 
on the number of bonds is predicted to be either linear, like a square root or 
logarithmic. 



PACS: 87.15 By, 82.20 Mj 



Introduction. Single molecule force spectroscopy has made it possible to measure the 
binding strength of a pair of receptor-ligand ("lock-key") molecules using vesicles [Q, the 
atomic force apparatus @-§J, or optical tweezers f5[ as transducers. Thus, the essential 
constituents mediating biological adhesion have become accessible to quantitative physical 
experiments ||. This experimental progress has fostered theoretical studies of the rupture 
of such pairs under dynamic loading. Thermal activation being a main contributing factor, 
Kramers-like descriptions of the rupture process with time-dependent potentials show that 
the rupture strength of such bonds depends on the loading rate |7|-§|. Such behavior has 
been found experimentally indeed 0,0. While unspecific theoretical models of the rupture 
process reveal generic features, molecular dynamic studies can address the details of the 
dynamics of the rupture of specific pairs fl,p~2]j. 

Adhesive contact and the rupture thereof often involves not just one but several molecular 
pairs of the same or different species [T3J. The equilibrium properties of the cooperative 
effects of such specific interactions are well studied both in theory 0-4T6| and in experiments 
|T7]Hl9| • Concerning the dynamics of rupture of such a contact under loading, detailed models 



for specific problems such as the peeling of a membrane |20"| , |2Tf or the rolling of leucocytes 
in shear flow [22] have been solved numerically to extract a critical tension or shear rate for 
rupture. However, it is inherently difficult to separate generic dynamical properties from 
specific ones using such intricate models. As an example for a generic property consider 
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the following question: How does the time and force necessary to break an adhesive contact 
under dynamic loading depend on the number of bonds initially present? 

The present study addresses this question within a simple model that extends work on 
the dynamic failure of a single bond to that of a whole patch involving several bonds of 
the same type. Quite generally, two different limiting cases must be distinguished. If the 
load is primarily concentrated on one bond at a time with relaxation of the load when the 
first bonds fails and subsequent loading of the next one, the rupture process basically is a 
sequence of similar single molecule events. The N bonds initially present then act in series. 
The rupture time will be ~ N whereas the force will exhibit a saw tooth-like pattern with 
a peak given by the rupture strength of a single bond. Such a behavior has been found and 
modeled in the related case of unfolding of proteins with several identical domains like titin 
23| . The main purpose of the present paper is to analyze the other case where the load is 



distributed (almost) uniformly among several bonds such that these bonds act in parallel. 
As more and more bonds rupture, the force on the remaining ones increases. This simple 
type of cooperativity leads to different scaling regimes for the rupture time and rupture 
force. 




FIG.l: Model geometry for the rupture of parallel bonds. Symbols are explained in the 
main text. 



Model. We model the rupture geometry generically as shown in Fig.l. One partner of 
the bond ("receptor") is confined to a substrate. The other ("ligand") is connected by a 
polymer to a transducer which is connected by an elastic element to a sled being pulled at 
velocity v. For simplicity, we model both the elasticity of the transducer and the polymers 
as Hookean springs with zero rest length and spring constants K t and K p , respectively. As 
long as a bond is intact the corresponding polymer is stretched to an extension x p which 
we assume to be the same for all intact bonds. The elongation of the transducer from its 
resting position is x t . Force balance on the transducer becomes NK p x p = K t x t where N 
is the number of intact bonds. Geometry dictates the time dependence x p (t) + x t (t) = vt. 
From these two relations, we find the time- dependent force on an intact bond as 



Following Bell ||13|| , we assume that the main effect of such a force is to introduce an instan- 



taneous, time-dependent dissociation rate k (t) according to 
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ko(t) = k exp[F b (t)xb/k B T], 



(2) 



where k is the dissociation rate in the absence of a force. The quantity x b is of the order of 
the distance between the minimum of the binding potential and the barrier and k B T is the 
product of Boltzmann's constant and temperature. 

We are mainly interested in the case of a soft transducer defined as K t <, K p . In this 
case, eq. (P shows that the force on a bond is inversely proportional to the number of intact 
bonds for all N(t). Hence, when a bond ruptures, the force on the remaining ones increases. 
We now discuss two different cases, irreversible and reversible bonds. In the former case, a 
bond, once ruptured, cannot rebind. Reversible bonds have a non-zero rebinding rate. 

Irreversible bonds. Initially N(t = 0) = N bonds are present. The rate equation for 
their time-dependent decrease is 

d t N = -N(t)k exp[F b (t)x b /k B T}. (3) 

We scale time with the dissociation rate in equilibrium fc according to r = tk Q . The rate 
equation in the case of a soft transducer then becomes 

8 T N= -Nexp[(jir/N] (4) 

with the loading parameter 

\l = K t x b v/k B Tk . (5) 

This simple rate equation seems not to have an analytical solution. However, its scaling 
behavior can be extracted by the following analysis. With the substitution u(r) = t/N one 
obtains 



d T u = u(1/t + 1) + u(exp[nu] — 1) (6) 

For small r, u(t) ~ t/Nq and the second term in (Qj can therefore be neglected. The 
solution U\{t) of the corresponding equation becomes u\{t) = re T /N and hence a purely 
exponential decay for the number of intact bonds, N(r) = N e~ T . This approximation 
breaks down for r >, T\ with T\ implicitly defined by 

max(l/ri, 1) = exp[^Wx(n)] — 1. (7) 

For r > n, we can then ignore both the first term and the in the second term of (^j). 
The corresponding equation d T u = wexpf/iw] is solved by 

E(jiu\) - E(fiu) =r-n (8) 

where E(x) = dx'e~ x ' /x' is the exponential integral and U\ = Wi(ti) is the cross over 
value of the first solution at the matching point t%. Hence the time necessary for complete 
rupture, r*, can be estimated by setting u(r*) = r*/N = oo which leads to 

T* = Ti+T 2 = T 1 + E( f m 1 ). (9) 
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Based on this approximative solution of @, three sub- regimes can be identified: 

(i) n < 1: In this case, the exponential decay holds till N(t) ~ 1. Physically, the rupture 
is then complete. In this trivial regime, where the loading is too small to affect the rupture 
process at all, the time required for rupture is 

r*~lniVo- (10) 

Note that the same result could have been obtained by analyzing the mean time required 
for the irreversible decay of N independent bonds under no force. 

(ii) 1 <, // <, N : In this regime, the exponential decay persists till T\ ~ \n(N / ft). At this 
time the number of bonds has reached N(ti) ~ fi. The remaining bonds decay according 
to (|5|) which leads to an additional time Ti of order 1 which is small compared to T\. Hence 
the whole rupture time in this regime is of order 

T* ~ HNof/i). (11) 

(iii) N < fi. In this case, the exponential decay applies till T\ ~ (N /fi) ln(/x/iVo). 
According to (§) the remaining time r 2 ~ (N / n)/ Iii(h/Nq) is smaller than T\. Hence the 
total rupture time is 

T*~(N /(j,)hy(ji/N ). (12) 

Thus we find for small loading rates that the rupture time is logarithmic in the number of 
bonds initially present whereas for large loading, this time becomes linear in N . For fixed N 
and increasing /i, the rupture time first is independent of fi. It then decays logarithmically 
in \i and finally becomes inversely proportional to fi. 

The force measured by the transducer is given by 

F t = N(t)K p x p (t) w (k B T/x b )i2T. (13) 

Thus, the total force experienced by the soft transducer is independent of the number of 
intact bonds and increases linearly in time. The dimensionless rupture force /* = fir* is 
thus given by 

/* ~/xlnJVo for fi < 1, 

~ filn(N Q //j,) for l<^<iVo, 

~ jVolnOx/iVo) for N < fi. (14) 

in the three regimes, respectively. 

Reversible bonds. So far, we have neglected the possibility that broken bonds can 
reform. Hence, rupture from a genuine equilibrium situation where bonds form, break, and 
rebind requires a refined description where we add a term for rebinding. We assume that one 
species of the receptor/ligand couple is limited to a total number Ni with N(t) molecules 
bound and N± — N(t) unbound whereas the other species is available in excess. The rate 
equation becomes 

d t N = -N(t)k exp\piT/N{t)] + kf(N x - N(t)), (15) 
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where we assume for simplicity that the rate kf for bond formation is not affected by the 
force. Without loading, the equilibrium number of bonds is 

N eq = 1 N 1 /(l+ 1 ) (16) 

where 7 = kf/k . As loading starts, the number of bonds decreases from this equilibrium 
value. With u(t) = r/N as before, we get 

d T u = u(1/t + 1 +7 - 7iVi«/V) + u(exp[/ju] - 1) (17) 

For ji = 0, this equation is solved by 

«o(r) = r/N eq} (18) 

which corresponds to the stationary equilibrium distribution. The loading term becomes 
relevant at a time r = T\ for which 

(exp|/xuo(ri)/JVi] - l)u ~ d T u = 1/N eq . (19) 

Two cases must then be distinguished: 

(i) For jj, < N eq , n ~ (Neg/ii) 1 / 2 . Up to this time, the loading has not significantly af- 
fected the number of bonds. The remaining time till all bonds are ruptured can be estimated 
to be of the same order as n using flTTp. Hence, 

r* ~ {N e J^\ (20) 

In this case, the rupture time increases as a square root of the equilibrium bonds present 
and decreases as a square root of the loading parameter. 

(ii) For N eq < fi, t\ ~ [N eq j \x) ln(fj,/N eq ), with a remaining time of the same order. Hence 
in this case, we recover the irreversible result ( |12"D with No replaced by N eq . 

Since in both cases the rupture time r* ~ ri, we get easily for the rupture force 

f* = fir* ~ (jiNeJ 1 ' 3 for fi<N eq} 

- N eq \n^/N eq ) for N eq < fi. (21) 

Stiff transducer. So far, we have considered the case of a soft transducer [K t < K p ) for 
which the force on a bond depends on the number of bonds. Another limiting case is a 
stiff transducer with K t >, N K p and K t > N eq K p for the case of irreversible and reversible 
rupture, respectively. According to eq. (|1|), the force on a bond then is (almost) independent 
of the number of bonds. Hence, the rupture time is only weakly dependent on the number 
of bonds. An analysis of the corresponding rate equations along similar lines as above shows 
for the irreversible case two subregimes with 

t* ~ In N for p, < 1, 

~ln/2//Z for ft > 1 (22) 

with a loading parameter 
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ji = K p x b v/k B Tk (23) 

dominated by the polymeric stiffness. For the dimensionless maximal force experienced by 
the transducer during the rupture process one finds 

/* ~ pN for p, < 1 , 

~ A" In Ji for ft > 1 (24) 

in the two cases. 

Similarly, for a stiff transducer and reversible bonds, one gets 

r*~ (InN^ 1 / 2 /^ 2 for ji < 1, 

~ In ///// for // > 1 (25) 

and for the dimensionless maximal force experienced by the transducer 

f* ~ fi 1/2 N eq for fl<l, 

~N eq \njj for jl>l. (26) 

Finally, there is a crossover regime for No j£q K p K t K p , where the pulling starts as in 
the soft case. As the number of intact bonds decreases towards the value A" = K t /K p , the 
denominator in ([J) becomes dominated by K t and the rupture process proceeds as for a stiff 
transducer. For the reversible case, it turns out that both the rupture time and the rupture 
force are dominated by the soft part. Hence the results p0|j21| ) apply for all K t < N eq K rj 
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For the irreversible case, analysis of the crossover regime is slightly more involved. The 
different scaling regimes for rupture time and force are shown in Fig. 2 without explicit 
derivation. 

Concluding perspective. Based on an analysis of rate equations, the comprehensive scaling 
analysis presented in this paper has revealed several different regimes for the rupture time 
and force of parallel molecular bonds under dynamic loading. The most distinctive regime 
is presumably the square root dependence of rupture time and force (|20|j2T| ) on loading rate 



and number of bonds derived for reversible bonds under small loading. Such a square root 
behavior on the loading rate is different from both the irreversible case and the dependence 
on loading rate for rupture of a single bond or bonds in series. An experimental result 
showing such an exponent could therefore be taken as a signature of breaking multiple 
parallel reversible bonds. Of course, it will be important to work with a model system 
where the number or density of bonds of at least one partner can be controlled in order to 
extract the dependence of rupture time and force on this crucial quantity. 

An obvious theoretical refinement of the present model would be to include fluctuations 
of the rupture time for individual bonds. Other ramifications can include allowing lateral 
interactions between the bonds, combining the simplistic Hookean transducer with a mem- 
brane patch with its own elasticity, or modeling the rupture process more delicately than 
done here to name just a few possibilities. It will be interesting to see how robust the 
scaling regimes derived in this paper will be under such modifications which can effectively 
lead to scenarios somewhere between the present "in parallel" case and the "in series" case 
described briefly in the introduction. Finally, it should be clear that in spite of - or rather 
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because of - the progress made in understanding the single bond behavior, the cooperative 
effects of several bonds under dynamic loading deserve further attention both in theory and 
in experiment. 

Acknowledgments: I thank E. Sackmann for a stimulating discussion and J. Shillcock for 
a critical reading of the manuscript. 
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FIG. 2: Dynamical phase diagram for (a) the dimensionless rupture time r* and (b) 
the dimensionless rupture force /* as a function of the two loading parameters fi (|5|) and 
/2 (^_f) in the case of irreversible rupture. In the region w, the rupture force is given by 
/• // lii( .\ 0/ ////). 
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